IFUM 650/FT-99 Bicocca-FT-99-37 
Out— of— equilibrium dynamics of large— 0^ QFT in finite volume 



C. Destri and E. Manfredini 

(a) Dipartimento di Fisica G. Occhialini, 
Universitd di Milano-Bicocca and INFN, sezione di Milano^'^ 
(b) Dipartimento di Fisica, Universitd di Milano 
and INFN, sezione di Milano^''^ 
(November 1999) 

Abstract 

The A(/)^ model in a finite volume is studied in the infinite limit both 
at equilibrium and out of equilibrium, with particular attention to certain 
fundamental features of the broken symmetry phase. The numerical solution 
of the dynamical evolution equations shows that the zero-mode quantum 
fluctuations cannot grow macroscopically large starting from microscopic ini- 
tial conditions. Thus we conclude that there is no evidence for a dynamical 
Bose-Einstein condensation, in the usual sense. On the other hand, out of 
equilibrium the long-wavelength fluctuations do scale with the linear size of 
the system, signalling dynamical infrared properties quite different from the 
equilibrium ones characteristic of the same approximation scheme. 
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I. INTRODUCTION 



In the last few years a great deal of attention has been paid to the study of interacting 
quantum fields out of equilibrium. There are, in fact, many interesting physical situations in 
which the standard S-matrix approach cannot give sensible information about the behavior 
of the system, because it evolves through a series of highly excited states (i.e., states of finite 
energy density). 

As an example consider any model of cosmological inflation: it is not possible to extract 
precise predictions on physical observables without including in the treatment the quantum 
back-reaction of the field on the space-time geometry and on itself 

On the side of particle physics, the ultra-relativistic heavy-ion collisions, scheduled in the 
forthcoming years at CERN-SPS, BNL-RHIC and CERN-LHC, are supposed to produce 
hadron matter at very high densities and temperatures; in such a regime the usual approach 
based on particle scattering cannot be considered a good interpretative tool at all. To 
extract sensible information from the theory new computational schemes are necessary, that 
go beyond the simple Feynman diagram expansion. The use of resummation schemes, like 
the Hartree-Fock (HF) approximation and the large limit (LN) 0, or the Hard 
Thermal Loop resummation for systems at finite temperature (HTL) can be considered 
a first step in this direction. They, in fact, enforce a sum over an infinite subset of Feynman 
diagrams that are dominant in a given region of the parameter space, where the simple 
truncation of the usual perturbative series at finite order cannot give sensible answers. 

Quite recently HF and LN have been used in order to clarify some dynamical aspects of 
the large (f)'^ theory. For the reader's benefit and to better motivate this work, we give a 
very short summary of the conclusions reached in previous works: i) the early time evolution 
is dominated by a so-called "linear regime" , during which the energy initially stored in one 
(or few) modes of the field is transferred to other modes via either parametric or spinodal 
unstabilities, resulting in a large particle production and a consequent dissipation for the 
initial condensate p|; ii) the linear regime stops at a time scale ti oc log(A~^) (where A is 
the quartic coupling constant), by which the effects of the quantum fluctuation become of 
the same order as the classical contribution and the dynamics turns completely non lin- 
ear and non perturbative iii) after the time ti the relaxation occurs via power laws 
with anomalous dynamical exponent [Q; iv) the asymptotic particle distribution, obtained 
as the result of the copious particle production at the expenses of the "classical" energy, 
is strongly non-thermal 1^,^]; and finally, v) at very large time scale, t ~ \/V (where V 
stands for the volume of the system), the non-perturbative and non-linear evolution might 
eventually produce the onset of a novel form of non-equilibrium Bose-Einstein condensation 
of the long-wavelength Goldstone bosons usually present in the broken symmetry phase of 
the model P,[lO|. Another very interesting result in concerns the dynamical Maxwell 
construction, which reproduces the fiat region of the effective potential in case of broken 
symmetry as asymptotic fixed points of the background evolution. Moreover, the LN ap- 
proximation scheme has been used to follow the evolution of a initial state characterized 
by an occupied spherical shell in momentum space (a spherical 'tsunami' |]11|), around a 



particular momentum fco , with the following results: i) in a theory where the symmetry 
is spontaneously broken at zero density, if we start with a finite density initial state with 
restored symmetry, the spinodal instabilities lead to a dynamical symmetry breaking; ii) 
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the evolution produces a re- arrangement of the particle distribution towards low momenta, 
signalling the onset of Bose condensation; iii) the equation of state of the asymptotic gas is 
ultra-relativistic (even if the distribution is not thermal) |Tl . 



In this article we present a detailed study, in finite volume, of dynamical evolution out 
of equilibrium for the scalar field in the large N limit. More precisely, we determine 
how such dynamics scales with the size of the periodic box containing the system in the 
case of uniform backgrounds. This is necessary to address questions like out-of-equilibrium 
symmetry breaking and dynamical Bose-Einstein condensation. 

In section |T| we define the model in finite volume, giving all the relevant notations and 
definitions. We also stress the convexity of the effective potential as an exact result, valid 
for the full renormalized theory in any volume. 

In section we derive the large approximation of the 0(iV)— invariant version of 
X((f)'^y model, according to the general rules of ref. [jl2|. In this derivation it appears evident 



the essential property of the N ^ oo limit of being a particular type of classical limit, so 
that it leads to a classical phase space, a classical hamiltonian with associated Hamilton's 
equations of motion [see eqs. ( |3.15| ), ( |3.16| ) and ( p. IT) )]. We then minimize the hamiltonian 
function(al) and determine the conditions for massless Goldstone bosons (i.e. transverse 
fluctuations of the field) to form a Bose-Einstein condensate, delocalizing the vacuum field 



expectation value (see also ref. |]T3|)- This necessarily requires that the width of the zero- 
mode fluctuations becomes macroscopically large, that is of the order of the volume. Only 
when the background takes one of the extremal values proper of symmetry breaking the 
width of the zero-mode fluctuations is of order L^/^, as typical of a free massless spectrum. 

The study of the lowest energy states of the model is needed for comparison with the 
results of the numerical simulations, which show that the zero-mode width ctq stays micro- 
scopic (that is such that ao/volume— > when the volume diverges) whenever it starts from 
initial conditions in which it is microscopic. Our results, in fact, show clearly the presence 
of a time scale tl, proportional to the linear size L of the system, at which finite volume 
effects start to manifest. We shall give a very simple physical interpretation of this time 
scale in section 1IID| . The important point is that after tl the zero mode amplitude starts 
decreasing, then enters an erratic evolution, but never grows macroscopically large. This re- 
sult is at odd with the interpretation of the linear late-time growth of the zero-mode width 
as a full dynamical Bose-Einstein condensation of Goldstone bosons, but is compatible with 
the "novel" form of BEC reported in P-[Tl . 

In fact we do find that the size of the low-lying widths at time tl is of order L, to be 
compared to the equilibrium situation where they would be of order L° in the massive case 
or of order L^/^ in the massless case. Perhaps the denomination "microscopic" should be 
reserved to this two possibilities. Therefore, since our initial condition are indeed microscopic 
in this restricted sense, we do observe in the out-of-equilibrium evolution a rapid transition 
to a different regime intermediate between the microscopic one and the macroscopic one 
characteristic of Bose-Einstein condensation. As we shall discuss more in detail later on, this 



fully agrees with the result found in \TU\, that the time-dependent field correlations vanish 
at large separations more slowly than for equilibrium free massless fields (as rather than 
r~^), but definitely faster than the equilibrium broken symmetry phase characterized by 
constant correlations at large distances. 

At any rate, when one considers microscopic initial conditions for the choice of bare 
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mass which corresponds to broken symmetry, the role itself of symmetry breaking is not 
very clear in the large N description of the out-of-equilibrium dynamics, making equally 
obscure the issues concerning the so-called quantum phase ordering |T^. This is because 
the limit — oo is completely saturated by gaussian states, which might signal the onset 
of symmetry breaking only developing macroscopically large fluctuations. Since such fluc- 
tuations do not appear to be there, the meaning itself of symmetry breaking, as something 
happening as times goes on and accompanied by some kind of phase ordering, is quite un- 
clear. We postpone to a companion work [|14| the discussion about the possibility of using 
more comprehensive approximation schemes, that include some non-gaussian features of the 
complete theory. As far as the large approximation is concerned, we underline that an 
important limitation of our approach, as well as of those of the references mentioned above, 
is in any case the assumption of a uniform background. Nonetheless, phenomena like the 
asymptotic vanishing of the effective mass and the dynamical Maxwell construction, taking 
place in this contest of a uniform background and large A^ expansion, are certainly very 
significant manifestations of symmetry breaking and in particular of the Goldstone theorem 
which applies when a continuous symmetry is broken. 

Finally, in section |^ we summarize the results presented in this article and we sketch 
some interesting open problems that we plan to study in forthcoming works. 



II. CUTOFF FIELD THEORY 

We consider the A^— component scalar field operator in a D— dimensional periodic box 
of size L and write its Fourier expansion as customary 

0(0;) =L-^/^ 5:0,6^^-, 4 = 0-. 

k 

with the wavevectors k naturally quantized: k = {27c /L)n, n G . The canonically conju- 
gated momentum tt has a similar expansion 

k 

with the commutation rules [4'k^'^-k'] — "^"^ifc*'^^"^- The introduction of a finite volume 
should be regarded as a regularization of the infrared properties of the model, which allows 
to "count" the different field modes and is needed especially in the case of broken symmetry. 
In fact, all the results we have summarized in section |I|, have been obtained simulating 
the system directly in infinite volume, where the evolution equations contain momentum 
integrals, that must be computed numerically by a proper, but nonetheless rather arbitrary, 
discretization in momentum space. Of course, the final result should be as insensitive as 
possible to the particular choice of the integration grid. In such a situation, the definition of 
a "zero" mode and the interpretation of its late time behavior might not be rigorous enough, 
unless, for some reason, it turns out that a particular mode requires a different treatment 
compared to the others. In order to understand this point, it is necessary to put the system 
in a finite volume (a box of size V); the periodic boundary conditions let us single out the 
zero mode in a rigorous way and thus we can carefully analyze its scaling properties with V 
and get some information on the infinite volume limit. 
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To regularize also the ultraviolet behavior, we restrict the sums over wavevectors to the 
points lying within the D— dimensional sphere of radius A, that is /c^ < A^, with TV = KL/2tt 
some large integer. Clearly we have reduced the original field-theoretical problem to a 
quantum-mechanical framework with finitely many (of order M^~^) degrees of freedom. 

The d)^ Hamiltonian reads 



+ 



H=\jd''x [tt^ + {dct^f + ml4>' + K{<t>y] = ^ E h - ^-fc + + ^b) <t>u- 

k 

+ ^ E i.(t>k{ 4>k2)i.'t>k,- (t>k,) ^i^,lk2+k,+k,fl (2-1) 
^1,^2,^3,^4 

where and Ab should depend on the UV cutoff A in such a way to guarantee a finite 



limit A — »• oo for all observable quantities. As is known |§|,|T5|, this implies triviality (that 
is vanishing of renormalized vertex functions with more than two external lines) for D > ?> 
and very likely also for D = ?>. In the latter case triviality is manifest in the one-loop 
approximation and in large— limit due to the Landau pole. For this reason we shall keep 
A finite and regard the model as an effective low-energy theory (here low-energy means 
practically all energies below Planck's scale, due to the large value of the Landau pole for 
renormalized coupling constants of order one or less). 



We shall work in the wavefunction representation where {ip\ \E') = ^^(v?) and 

d 



d 



while for > (in lexicographic sense) 

1 1 f d d \ 

(0±fc^)(c^) = -= (cpfc ± ? ip_,) ^{cp) , {7T±k^){cp) = -= -t—- ± ^(C^) 



Notice that by construction the variables ipi^ are all real. Of course, when either one of 
the cutoffs are removed, the wave function \l/(c^) acquires infinitely many arguments and 
becomes what is usually called a wavefunctional. 

In practice, the problem of studying the dynamics of the 0^ field out of equilibrium 
consists now in trying to solve the time-dependent Schroedinger equation given an initial 
wavefunction \I^(</',t = 0) that describes a state of the field far away from the vacuum. 
By this we mean a non-stationary state that, in the infinite volume limit L ^ oo, would 
lay outside the particle Fock space constructed upon the vacuum. This approach could be 
generalized in a straightforward way to mixtures described by density matrices, as done, for 



instance, in |16-18|. Here we shall restrict to pure states, for sake of simplicity and because 
all relevant aspects of the problem are already present in this case. 

It is by now well known |^ that perturbation theory is not suitable for the purpose 
stated above. Due to parametric resonances and/or spinodal instabilities there are modes 
of the field that grow exponentially in time until they produce non-perturbative effects 
for any coupling constant, no matter how small. On the other hand, only few, by now 
standard, approximate non-perturbative schemes are available for the 0^ theory, and to 
these we have to resort after all. We shall consider here only the large N expansion to 
leading order, remanding to another work the definition of a time-dependent Hartree-Fock 
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(tdHF) approach (a generalization of the treatment given, for instance, in [^]). In fact 
these two methods are very closely related, as shown in where several techniques to 



derive reasonable dynamical evolution equations for non-equilibrium 0^ are compared. 

We close this section by stressing that the introduction of both a UV and IR cutoff 
allows to easily derive the well-known rigorous result concerning the flatness of the effective 
potential. In fact Vcs{<f)) is a convex analytic function in a finite neighborhood of = 0, as 
long as the cutoffs are present, due to the uniqueness of the ground state. In the infrared 
limit L — »• oo, however, Kfr(0) might flatten around = 0. Of course this possibility would 
apply in case of spontaneous symmetry breaking, that is for a double-well classical potential. 
This is a subtle and important point that will play a crucial role later on, even if the effective 
potential is relevant for the static properties of the model rather than the dynamical evolution 
out of equilibrium that interests us here. In fact, the dynamical evolution in QFT is governed 
by the OTP effective action [^,0 and one might expect that, although non-local in time, 
it asymptotically reduces to a multiple of the effective potential for trajectories of (pit) with 
a fixed point at infinite time. In such case there should exist a one-to-one correspondence 
between fixed points and minima of the effective potential. 



III. LARGE EXPANSION AT LEADING ORDER 

A. Definitions 

In this section we consider the standard non-perturbative approach to the 0*^ model 
which is applicable also out of equilibrium, namely the large method as presented in 



1 2^ . However we shall follow a different derivation which makes the gaussian nature of the 
N oo limit more explicit. 

It is known that the theory described by the Hamiltonian ( |2.1| ) is well behaved for large 
N, provided that the quartic coupling constant Ab is rescaled with For example, it 

is possible to define a perturbation theory, based on the small expansion parameter 
in the framework of which one can compute any quantity at any chosen order in 
From the diagrammatic point of view, this procedure corresponds to a resummation of the 
usual perturbative series that automatically collects all the graphs of a given order in 
together Moreover, it has been established since the early 80's that the leading order 



approximation (that is the strict limit oo) is actually a classical limit |12], in the sense 
that there exists a classical system (i.e., a classical phase space, a Poisson bracket and a 
classical Hamiltonian) whose dynamics controls the evolution of all fundamental quantum 
observables, such as field correlation functions, in the iV — >■ cx) limit. For instance, from the 
absolute minimum of the classical Hamiltonian one reads the energy of the ground state, 
while the spectrum is given by the frequencies of small oscillations about this minimum, etc. 
etc.. We are here interested in finding an efficient and rapid way to compute the quantum 
evolution equations for some observables in the N oo limit, and we will see that this task 
is easily accomplished just by deriving the canonical Hamilton equations from the large 
classical Hamiltonian. 



Following Yaffe [|12| , we write the quantum mechanical hamiltonian as 

H = Nh{A,C) (3.1) 
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in terms of the square matrices A, C with operator entries (-07^ is the canonical momentum 
conjugated to the real mode (pf,) 



A-kk' — jyV'fc' Vk' ; Ckk' = jy'^fc' ''^k' 

These are example of "classical" operators, whose two-point correlation functions factorize 
in the N ^ 00 limit. This can be shown by considering the coherent states 



'^z,q,p{^) = C{z) exp 



2N 



Y^^kk'i^k- Qk) ■ i^Pk'-Qk' 



kk' 



(3.2) 



where the complex symmetric matrix z has a positive definite real part while and q^. are 
real and coincide, respectively, with the coherent state expectation values of v^k and (pf.. As 
these parameters take all their possible values, the coherent states form an overcomplete set 
in the cutoff Hilbert space of the model. The crucial property which ensures factorization 
is that they become all orthogonal in the N ^ 00 limit. Moreover one can show that 
the coherent states parameters form a classical phase space with Poisson brackets 



{Qk,Pk'} 



P.B. 



5kk'5'^ , {Wkk' ,Vgg'} 



P.B. 



SkqSk'g' + Skq'S, 



kq'Ok'q 



where w and v reparametrize z as z 



+ iv. It is understood that the dimensionality 



of the vectors q^^ and is arbitrary but finite [that is, only a finite number, say n, of pairs 
{ip\ , Tu\) may take a nonvanishing expectation value as 00]. 

Once applied to the classical operators Akk' and Ckk' the large factorization allow to 
obtain the classical hamiltonian by simply replacing A and C in eq. (|3 . 1|) by the coherent 
expectation values 



{Akk') — Qk ■ Qk' + '^kk' , {Ckk') 



Pk-Pk' + [vwvjkk' + -{w- 



Ikk' 



In our situation, having assumed a uniform background expectation value for 0, we have 
Qk — Pk — ^ for all A; 7^ 0; moreover, translation invariance implies that w and v are diagonal 
matrices, so that we may set 



Wkk' 



k '^kk' 1 Vkk' 



Sk r 

Okk' 

0"fc 



in term of the canonical couples (cTfc, Sk) which satisfy {ak , Sfc'jp.B. = ^kk'- Notice that the 
cTfc are just the widths (rescaled by A^~^/^) of the 0{N) symmetric and translation invariant 
gaussian coherent states. 

Thus we find the classical hamiltonian 

2 



h 



cl 



2 2\ 

Qo) 



2r 



+ 



4L^ 



where by Hamilton's equations of motion Pq = and Sk = crk- The corresponding conserved 
energy density S = L~^hc\ may be written 



8 = T + V 



r 



li^ 1 
2^ +2l^ 



E^^ 



V 



2L^ 



1 



+ + S) 



1 



(3.3) 



E-, 

k 
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where = L'^/'^q^ and V is the 0(A^)— invariant quartic potential regarded as a function 
of </)^, that is V{z) = ^m^z + \XhZ^- It is worth noticing that eq. ( P73| ) would apply as is to 
generic V{z). 



B. Static properties 



Let us consider first the static aspects embodied in the effective potential ^^(0), that 
is the minimum of the potential energy V at fixed (f). We first define in a precise way the 
unbroken symmetry phase, in this large context, as the case when Kff(0) has a unique 
minimum at = in the limit of infinite volume. Minimizing V w.r.t. (Tk yields 



2VW+JP 



2 I , 



ml + 2V'{cl> 
ml + X^cj)'^ + 



(3.4) 



that is the widths characteristic of a free theory with self-consistent mass M fixed by the gap 
equation. By the assumption of unbroken symmetry, when (j) = and at infinite volume M 
coincides with the equilibrium mass m of the theory, that may be regarded as independent 
scale parameter. Since in the limit L ^ oo sums are replaced by integrals 



we obtain the standard bare mass parameterization 



m. 



- Ab/D("^^ A) , loiz, A)= f 

Jk 



and the renormalized gap equation 

= + A 0^ + A [/D(M^ A) - A) 
which implies, when D = 3, 



finite 



Ab = A 1 



A , 2A 
— log 



,-1 



m\/e 



(3.5) 



(3.6) 



(3.7) 



with a suitable choice of the finite part. No coupling constant renormalization occurs instead 
when D = 1. The renormalized gap equation ( |3.6|) may also be written quite concisely 



M2 



m 



A(M) A(m 
in terms of the one-loop running couplings constant 

A 



+ 0^ 



(3.8) 



A(/i) = A 



1 ^ 

log — 

m 



A(m) = A, A(2Ae-^/2) = Ab 
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It is the Landau pole in A(2Ae~^/^) that actually forbids the limit A — oo. Hence we must 
keep the cutoff finite and smaller than Apoie, so that the theory does retain a slight inverse- 
power dependence on it. At any rate, there exists a very wide window where this dependence 
is indeed very weak for couplings of order one or less, since Apoie = |mexp(l/2+87r^/A) 3> m. 
Moreover, we see from eq. (p.8|) that for V^|0| much smaller than the Landau pole there 
are two solutions for M, one "physical", always larger than m and of the same order of 
m + \/X\(j)\, and one "unphysical" , close to the Landau pole. 

One can now easily verify that the effective potential has indeed a unique minimum in 
= 0, as required. In fact, if we assign arbitrary 0— dependent values to the widths 0"^, 
(minus) the effective force reads 

|,y»,{..W})^M= ^-.1:1110 (3.9) 

and reduces to 0* when the widths are extremal as in eq. (|3.4|) ; but is positive for 
unbroken symmetry and so = is the unique minimum. 

We define the symmetry as broken whenever the infinite volume has more than one 
minimum. Of course, as long as L is finite, V^g has a unique minimum in = 0, because of 
the uniqueness of the ground state in Quantum Mechanics, as already discussed in section 
Let us therefore proceed more formally and take the limit L oo directly on the potential 
energy V. It reads 

where we write for convenience the tree-level potential V in the positive definite form V{z) = 
\\h{z + rrib/Ab)^. V is now the sum of two positive definite terms. Suppose there exists a 

— 2 

configuration such that V{4> + S) = and the first term in V is at its minimum. Then this 
is certainly the absolute minimum of V. This configuration indeed exists at infinite volume 
when D = 3: 

""'^2^' = ^b = -Ab[t^' + /3(0,A)] (3.10) 

where the nonnegative v should be regarded as an independent parameter fixing the scale 
of the symmetry breaking. It replaces the mass parameter m of the unbroken symmetry 
case: now the theory is massless in accordance with Goldstone theorem. On the contrary, 
if D = 1 this configuration is not allowed due to the infrared divergences caused by the 
massless nature of the width spectrum. This is just the standard manifestation of Mermin- 
Wagner-Coleman theorem that forbids continuous symmetry breaking in a two-dimensional 



space-time 23 



At finite volumes we cannot minimize the first term in V since this requires (Jq to diverge, 

— 2 

making it impossible to keep V{(f) + S) = 0. In fact we know that the uniqueness of the 



ground state with finitely many degrees of freedom implies the minimization equations ( p.4|) 
to hold always true with a strictly positive. Therefore, broken symmetry should manifest 
itself as the situation in which the equilibrium value of is a positive definite function of 
L which vanishes in the L oo limit. 
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We can confirm this qualitative conclusion as follows. We assume that the bare mass has 
the form given in eq. ( p.lOD and that (j) = v"^ too. Minimizing the potential energy leads 
always to the massive spectrum, eq. (p.4|), with the gap equation 



^ ~ 2I?M ^ 21? ~^/WTW ~ 8^ ^ ■ ^ 

If > does not vanish too fast for large volumes, or stays even finite, then the sum on 
the modes has a behavior similar to the corresponding infinite volume integral: there is a 
quadratic divergence that cancels the infinite volume contribution, and a logarithmic one 
that renormalizes the bare coupling. The direct computation of the integral would produce 
a term containing the M^log(A/M). This can be split into M^[log(A/f) — log(M/f)] by 
using V as mass scale. The first term renormalizes the coupling correctly, while the second 
one vanishes if vanishes in the infinite volume limit. 
When L oo, the asymptotic solution of ( |3.11|) reads 



M= L-' + h.o.t. 

that indeed vanishes in the limit. Note also that the exponent is consistent with the as- 
sumption made above that M vanishes slowly enough to approximate the sum over k ^ 
with an integral with the same M. 

Let us now consider a state whose field expectation value (f) is different from f^. If 
4> > V, the minimization equations ( |3.4|) leads to a positive squared mass spectrum for the 
fluctuations, with given self-consistently by the gap equation. On the contrary, as soon 
as < w^, one immediately see that a positive cannot solve the gap equation 

if we insist on the requirement that ctq not be macroscopic. In fact, the r.h.s. of the previous 
equation is negative, no matter which positive value for the effective mass we choose, at 
least for L large enough. But nothing prevent us to consider a static configuration for which 
the amplitude of the zero mode is macroscopically large (i.e. it rescales with the volume 
L^). Actually, if we choose 

4 = "=-*=+ ' 



L3 ^ 2L3M 

we obtain the same equation as we did before and the same value for the potential, that is 
the minimum, in the limit L —>■ oo. Note that at this level the effective mass M needs not 
to have the same behavior in the L ^ oo limit, but it is free of rescaling with a different 
power of L. We can be even more precise: we isolate the part of the potential that refers to 
the zero mode width ctq (S' does not contain the o"o contribution) 
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and we minimize it, keeping 4^ fixed. The minimum is attained at t = ctq/L^ solution of 
tlie cubic equation 

Abt^ + aAbt^ - \L-^ = 

where a = (f)^ — v'^ + H' — (0, A). Note that Aba depends on L and it has a finite hmit in 
infinite volume: X{(f)^ — v'^). The solution of the cubic equation is 

Abt = X^{v' - 4>^) + \[L\v^ - + h.o.t. 

from which the effective mass can be identified as proportional to . The stability equa- 
tions for all the other modes can now be solved by a massive spectrum, in a much similar 
way as before. 

Since cxo is now macroscopically large, the infinite volume limit of the at distribution 
(that gives a measure of the transverse fluctuations in the 0{N) model) develop a 5— like 
singularity, signalling a Bose condensation of the Goldstone bosons: 

al = iv' - 6^^\k) + ^ (3.12) 



2 

At the same time it is evident that the minimal potential energy is the same as when cp = v 
that is the effective potential flattens, in accord with the Maxwell construction. 

Eq. ( p. 121) corresponds in configuration space to the 2— point correlation function 



,2 



{4>{x) ■ 4>{y)) _ 3.2 , f d'^k .2_it.„-,) 



where A£)[x — y) is the massless free-field equal-time correlator, while 



(3.13) 



C(^') = Q{v^ - ^') + 6(0' - v^) = max(t;2, 0^) (3.14) 

— 2 — 2 

This expression can be extended to unbroken symmetry by setting in that case C{4> ) = (f) . 



Quite evidently, when eq. (|3.14|) holds, symmetry breaking can be inferred from the limit 
\x — y\ oo, if clustering is assumed pi| , p5[ , since /S.£,{x — y) vanishes for large separations. 



Of course this contradicts the infinite volume limit of the finite-volume definition, cf) = 
liuiN^oo N~^/'^ {4){x)) , except at the extremal points (j) = v'^. 

In fact the L oo limit of the finite volume states with 0^ < v'^ do violate cluster- 
ing, because they are linear superpositions of vectors belonging to superselected sectors and 
therefore they are indistinguishable from statistical mixtures. We can give the following in- 
tuitive picture for large A^. Consider any one of the superselected sectors based on a physical 
vacuum with (f) = v"^. By condensing a macroscopic number of transverse Goldstone bosons 
at zero-momentum, one can build coherent states with rotated (f>. By incoherently averaging 
over such rotated states one obtains new states with field expectation values shorter than v 
by any prefixed amount. In the large A^ approximation this averaging is necessarily uniform 
and is forced upon us by the residual 0{N — 1) symmetry. 
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C. Out of equilibrium dynamics 



We now turn to the dynamics out of equilibrium in this large context. It is governed 
by the equations of motion derived from the total energy density £ in eq. (|3.3|), that is 



^4 



where the generally time-dependent effective squared mass is given by 



+ At 



0' + S-/D(m^A) 



in case of unbroken symmetry and 



(3.15) 



(3.16) 



(3.17) 



for broken symmetry in D = 3. 

At time zero, the specific choice of initial conditions for ak that give the smallest energy 
contribution, that is 



cTfc = , ai. 



1 



(3.18) 



turns eq. ( ^.161 ) into the usual gap equation (p.4[). For any value of (f) this equation has one 
solution smoothly connected to the value M = m at = 0. Of course other initial conditions 
are possible. The only requirement is that the corresponding energy must differ from that 
of the ground state by an ultraviolet finite amount, as it occurs for the choice (|3.18|) . In fact 
this is guaranteed by the gap equation itself, as evident from eq. (^.9|): when the widths ak 
are extremal the effective force is finite, and therefore so are all potential energy differences. 
This simple argument needs a refinement in two respects. 

Firstly, in case of symmetry breaking the formal energy minimization w.r.t. ak leads 
always to eqs. ( p.l^ ), but these are acceptable initial conditions only if the gap equation 
that follows from eq. ( |3.17 ) in the L —>■ oo limit, namely 



0^ 



(3.19) 



admits a nonnegative, physical solution for M^. 

Secondly, ultraviolet finiteness only requires that the sum over k in eq. (|3.9|) be finite 



and this follows if eq. (|3.18| ) holds at least for k large enough, solving the issue raised in the 
first point: negative are allowed by imposing a new form of gap equation 



-,2 



E 



1 



(3.20) 



where all ak with k'^ < |M^| are kept free (but all by hypothesis microscopic) initial condi- 
tions. Of course there is no energy minimization in this case. To determine when this new 
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form is required, we observe that, neglecting the inverse-power corrections in the UV cutoff 
we may write eq. (|3.19 ) in the following form 



J^ = 4>'-v' (3.21) 
A(M) ^ ^ 



There exists a positive solution smoothly connected to the ground state, = f ^ and 
= 0, only provided (f) > v"^. So, in the large N limit, as soon as we start with cf) < v"^, 
we cannot satisfy the gap equation with a positive value of M^. 

Once a definite choice of initial conditions is made, the system of differential equations 
( p.l5| ), (|3.16|) or ( p. 17] ) can be solved numerically with standard integration algorithms. This 



has been already done by several authors P-10|, working directly in infinite volume, with 



the following general results. In the case of unbroken symmetry it has been established that 
the ak corresponding to wavevectors k in the so-called forbidden bands with parametric 
resonances grow exponentially in time until their growth is shut off by the back-reaction. 
For broken symmetry it is the region in fc— space with the spinodal instabilities caused by 
an initially negative M^, whose widths grow exponentially before the back-reaction shutoff. 
After the shutoff time the effective mass tends to a positive constant for unbroken symmetry 
and to zero for broken symmetry (in D=3), so that the only width with a chance to keep 
growing indefinitely is ctq for broken symmetry. 

Of course, in all these approaches the integration over modes in the back-reaction E 
cannot be done exactly and is always replaced by a discrete sum of a certain type, depending 
on the details of the algorithms. Hence there exists always an effective infrared cutoff, albeit 
too small to be detectable in the numerical outputs. A possible troublesome aspect of this 
is the proper identification of the zero-mode width (Tq. Even if a (rather arbitrary) choice 
of discretization is made where a ctq appears, it is not really possible to determine whether 
during the exponential growth or after such width becomes of the order of the volume. Our 
aim is just to answer this question and therefore we perform our numerical evolution in finite 
volumes of several growing sizes. Remanding to the appendix for the details of our method, 
we summarize our results in the next subsection. 



D. Numerical results 

After a careful study in D = 3 of the scaling behavior of the dynamics with respect 
to different values of L, the linear size of the system, we reached the following conclusion: 
there exist a L— dependent time, that we denote by tl, that splits the evolution in two parts; 
for t < Tl, the behavior of the system does not differ appreciably from its counterpart at 
infinite volume, while finite volume effects abruptly alter the evolution as soon as t exceeds 
Tl] moreover 

• Tl is proportional to the linear size of the box L and so it rescales as the cubic root of 
the volume. 

• Tl does not depend on the value of the quartic coupling constant A, at least in a first 
approximation. 
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The figures show the behavior of the width of the zero mode ctq (see Fig. [^), of the 
squared effective mass (see Fig. |^ ) and of the back-reaction E (see Fig. ||), in the 
more interesting case of broken symmetry. The initial conditions are chosen in several 
different ways (see the appendix for details), but correspond to a negative at early times 
with the initial widths all microscopic, that is at most of order L^^^. This is particularly 
relevant for the zero-mode width ctq, which is instead macroscopic in the lowest energy 
state when (j) < t>^, as discussed above. As for the background, the figures are relative to 
the simplest case (f> = = (f), but we have considered also initial conditions with > 0, 
reproducing the "dynamical Maxwell construction" observed in ref. [|^. At any rate, for the 
purposes of this work, above all it is important to observe that, due to the quantum back- 
reaction, rapidly becomes positive, within the so-called spinodal time [|-10], and then. 



for times before tl, the weakly dissipative regime takes place where oscillates around zero 
with amplitude decreasing as and a frequency fixed by the largest spinodal wavevector, 
in complete agreement with the infinite-volume results |jlO|- Correspondingly, after the 
exponential growth until the spinodal time, the width of the zero-mode grows on average 
linearly with time, reaching a maximum for t ^ r^. Precisely, (Tq performs small amplitude 
oscillations with the same frequency of around a linear function of the form A + Bt, 
where A, B ^ X~^^'^ (see Fig. ^), confirming what already found in refs. then quite 

suddenly it turns down and enters long irregular Poincare-like cycles. Since the spinodal 
oscillation frequency does not depend appreciably on L, the curves of ctq at different values 
of L are practically identical for t < t^. After a certain number of complete oscillations, 
a number that scales linearly with L, a small change in the behavior of (see Fig. ^ 
determines an inversion in ctq (see Fig. ^), evidently because of a phase crossover between 
the two oscillation patterns. Shortly after tl dissipation practically stops as the oscillations 
of stop decreasing in amplitude and become more and more irregular, reflecting the 
same irregularity in the evolution of the widths. 

We can give a straightforward physical interpretation for the presence of the time scale 
Tl- As shown in [rU|, long after the spinodal time ti, the effective mass oscillates around 
zero with a decreasing amplitude and affects the quantum fluctuations in such a way that 
the equal-time two-point correlation function contains a time-dependent non-perturbative 
disturbance growing at twice the speed of light. This is interpreted in terms of large numbers 
of Goldstone bosons equally produced at any point in space (due to translation invariance) 
and radially propagating at the speed of light. This picture applies also at finite volumes, 
in the bulk, for volumes large enough. Hence, due to our periodic boundary conditions, 
after a time exactly equal to L/2 the forward wave front meets the backward wave front at 
the opposite point with respect to the source, and the propagating wave starts interfering 
with itself and heavily changes the dynamics with respect to that in infinite volume. This 
argument leads us to give the value of tt for the proportionality coefficient between and 
L/27r, prevision very well verified by the numerical results, as can be inferred by a look at 
the figures. 

The main consequence of this scenario is that the linear growth of the zero-mode width 
at infinite volume should not be interpreted as a standard form of Bose-Einstein Conden- 
sation (BEC), occurring with time, but should be consistently considered as "novel" form 
of dynamical BEC, as found by the authors of [|1^]. In fact, if a macroscopic condensation 
were really there, the zero mode would develop a 6 function in infinite volume, that would 
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be announced by a width of the zero mode growing to values 0(L^/^) at any given size 
L. Now, while it is surely true that when we push L to infinity, also the time tl tends to 
infinity, allowing the zero mode to grow indefinitely, it is also true that, at any fixed though 
arbitrarily large volume, the zero mode never reaches a width 0(L'^/^), just because tl oc L. 
In other words, if we start from initial conditions where ctq is microscopic, then it never 
becomes macroscopic later on. 

On the other hand, looking at the behavior of the mode functions of momenta k = 
(27r/L)n for n fixed but for different values of L, one realizes that they obey a scaling 
similar to that observed for the zero-mode: they oscillate in time with an amplitude and a 
period that are 0{L) (see fig. |^ and §). Thus, each mode shows a behavior that is exactly 
half a way between a macroscopic amplitude [i.e. 0(L^/^)] and a usual microscopic one [i.e. 
at most 0(L^/^)]. This means that the spectrum of the quantum fiuctuations at times of 
the order of the diverging volume can be interpreted as a massless spectrum of interacting 
Goldstone modes, because their power spectrum develops in the limit a l//c^ singularity, 
rather than the 1/k pole typical of free massless modes. As a consequence the equal-time 
field correlation function [see eq. ( |3.13| )] will fall off as \x — y\^^ for large separations smaller 
only than the diverging elapsed time. This is in accord with what found in [0, where the 
same conclusion where reached after a study of the correlation function for the scalar field 
in infinite volume. 

The fact that each mode never becomes macroscopic, if it started microscopic, might be 
regarded as a manifestation of unitarity in the large approximation: an initial gaussian 
state with only microscopic widths satisfies clustering and clustering cannot be spoiled by 
a unitary time evolution. As a consequence, in the infinite-volume late-time dynamics, the 
zero-mode width do does not play any special role and only the behavior of as ^ is 
relevant. As already stated above, it turns out from our numerics as well as from refs. [p|-[TT|l 
that this behavior is of a novel type characteristic both of the out-of-equilibrium dynamics 
and of the equilibrium finite-temperature theory [^ , with oc 1/k. 

A final comment should be made about the periodic boundary conditions used for these 
simulations. This choice guarantees the translation invariance of the dynamics needed to 
consider a stable uniform background. If we had chosen other boundary conditions (Dirichlet 
or Neumann, for instance), the translation symmetry would have been broken and an uniform 
background would have become non-uniform pretty soon. Of course, we expect the bulk 
behavior to be independent of the particular choice for the boundary conditions in the 
infinite volume limit, even if a rigorous proof of this statement is still lacking. 



IV. DISCUSSION AND PERSPECTIVES 

In this work we have presented a rather detailed study of the dynamical evolution out of 
equilibrium, in finite volume (a cubic box of size L in 3D), for the 0^ QFT in the large A^ 
limit. For comparison, we have also analyzed some static characteristics of the theory both 
in unbroken and broken symmetry phases. 

We have reached the conclusion, based on strong numerical evidence, that the linear 
growth of the zero-mode quantum fluctuations, observed already in the large A^ approach of 
refs. [PHTTH, may be consistently interpreted as a "novel" form of dynamical Bose-Einstein 
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condensation, different from the traditional one in finite temperature field theory at equilib- 
rium. In fact, in finite volume, (Tq never grows to 0(L^/^) if it starts from a microscopic value, 
that is at most of order L^/^. On the other hand all long- wavelength fluctuations rapidly 
become of order L, signalling a novel infrared behavior quite different from free massless 
fields at equilibrium [recall that the large N approximation is of mean field type, with no 
direct interaction among particle excitations]. This is in agreement with the properties of 
the two-point function determined in [|T^]. 

The numerical evidence for the linear dependence of on L is very strong, and the qual- 
itative argument given in the previous section clearly explains the physics that determines 
it. Nonetheless a solid analytic understanding of the detailed (quantitative) mechanism 
that produces the inversion of do around tl and its subsequent irregular behavior, is, at 
least in our opinion, more difficult to obtain. One could use intuitive and generic argu- 
ments like the quantization of momentum in multiples of 2n/L, but the evolution equations 
do not have any simple scaling behavior towards a universal form, when mass dimensions 
are expressed in multiples of 27t/L and time in multiples of L. Moreover, the qualitative 
form of the evolution depends heavily on our choice of initial conditions. In fact, before 
finite volumes effects show up, the trajectories of the quantum modes are rather complex 
but regular enough, having a small-scale quasi-periodic almost mode-independent motion 
within a large-scale quasi-periodic mode-dependent envelope, with a very delicate resonant 
equilibrium [Cfr. Fig. 1 and 7]. Apparently [Cfr. Fig. 5 and 6], it is a sudden small beat 
that causes the turn around of the zero-mode and of the other low-lying modes (with many 
thousands of coupled modes, it is very difficult for the delicate resonant equilibrium to fully 
come back ever again), but we think that a deeper comprehension of the non-linear coupled 
dynamics is needed in order to venture into a true analytic explanation. 

On the other hand it is not difficult to understand why does not depend appreciably 
on the coupling constant: when finite- volume effects first come in, that is when the wave 
propagating at the speed of light first starts to interfere with itself, the quantum back- 
reaction AS has settled on values of order 1, because the time tl — L/2 is much greater 
than the spinodal time ti. The slope of the linear envelope of the zero mode does depend 
on A because it is fixed by the early exponential growth. Similarly, it is easy to realize that 
the numerical integrations of refs. [p|-p]] over continuum momenta correspond roughly to an 
effective volume much larger than ours, so that the calculated evolution remained far away 
from the onset of finite- volume effects. 

The main limitation of the large N approximation, as far as the evolution of the widths 
cTfc is concerned, is in its intrinsic gaussian nature. In fact, one might envisage a scenario in 
which, while gaussian fluctuations stay microscopic, non-gaussian fluctuations grow in time 
to a macroscopic size. Therefore, in order to clarify this point and go beyond the gaussian 
approximation, we are going to consider, in a forthcoming work |]14|, a time-dependent HF 
approximation capable in principle of describing the dynamics of non-gaussian fluctuation 
of a single scalar field with 0^ interaction. 

Another open question concerns the connection between the minima of the effective 
potential and the asymptotic values for the evolution of the background, within the simplest 



gaussian approximation. As already pointed out in |T0|, a dynamical Maxwell construction 
occurs for the 0{N) model in infinite volume and at leading order in in case of broken 
symmetry, in the sense that any value of the background within the spinodal region can 
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be obtained as large time limit of the evolution starting from suitable initial conditions. It 
would be very enlightening if we could prove this "experimental" result by first principles 
arguments, based on CTP formalism. Furthermore, preliminary numerical evidence |T^ 
suggests that something similar occurs also in the Hartree approximation for a single field, 
but a more thorough and detailed analysis is needed. 

It would be interesting also to study the dynamical realization of the Goldstone paradigm, 
namely the asymptotic vanishing of the effective mass in the broken symmetry phases, in 
different models; this issue needs further study in the 2D case where it is known that 
the Goldstone theorem is not valid. 
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APPENDIX A: DETAILS OF THE NUMERICAL ANALYSIS 

We present here the precise form of the evolution equations for the field background and 
the quantum mode widths, which control the out-of-equilibrium dynamics of the 0^ model 
in finite volume at the leading order in the 1/A^ approach, as described in sections [III C| . We 
restrict here our attention to the tridimensional case. 

Let us begin by noticing that each eigenvalue of the Laplacian operator in a ?)D finite 

volume is of the form A;^ = where n is a non-negative integer obtained as the sum 

of three squared integers, n = + + n^. Then we associate a degeneracy factor gn to 
each eigenvalue, representing the number of different ordered triples (n^., Uy, Hz) yielding the 
same n. One may verify that Qn takes on the continuum value of 47r/c^ in the infinite volume 

limit, where k = (^^^ n is kept fixed when L — oo. 

Now, the system of coupled ordinary differential equations is, in case of the large N 
approach. 



0, 



(Jr. 



4(7,3 







(Al) 



where the index n ranges from to A/"^, M = AL/27T and M^(t) is defined by the eq. 
( p.l6| ) in case of unbroken symmetry and by eq. ( p.l7| ) in case of broken symmetry. The 
back-reaction E reads, in the notations of this appendix 



1 

Id 



n=0 
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Technically it is simpler to treat an equivalent set of equations, which are formally linear 
and do not contain the singular Heisenberg term oc a~^. This is done by introducing the 
complex mode amplitudes = (T„exp(i^„), where the phases 9n satisfy o"^^„ = 1. Then we 
find a discrete version of the equations studied for instance in ref 



+ — n + M^ 



dt^ ' y^J - ' - ^-r. ^, ^ ^T.^9n\Zn\' (A2) 

subject to the Wronskian condition 

One realizes that the Heisenberg term in cr„ corresponds to the centrifugal potential for the 
motion in the complex plane of Zn- 

Let us now come back to the equations ([A^). To solve these evolution equations, we 
have to choose suitable initial conditions respecting the Wronskian condition. In case of 
unbroken symmetry, once we have fixed the value of and its first time derivative at initial 
time, the most natural way of fixing the initial conditions for the z„ is to require that they 
minimize the energy at t = 0. We can obviously fix the arbitrary phase in such a way to 
have a real initial value for the complex mode functions 

- ^(«' - 'vf 

where fi„ = \Jk'^ + M'^ifi). The initial squared effective mass M^(t = 0), has to be deter- 
mined self-consistently, by means of its definition ( |3.16|) . 

In case of broken symmetry, the gap equation is a viable mean for fixing the initial 
conditions only when lies outside the spinodal region [see eq ( |3.21| )]; otherwise, the gap 
equation does not admit a positive solution for the squared effective mass. In that case, 
we have to resort to other methods, in order to choose the initial conditions. Following the 
discussion presented in [III Q , one possible choice is to set a\ = ^— ^==== for k"^ < |M^| in 



eq. (|3.2CI|) and then solve the corresponding gap equation ( p. 201 ). An other acceptable choice 



would be to solve the gap equation (|3.2CI| ), once we have set a massless spectrum for all the 
spinodal modes but the zero mode, which is started from an arbitrary, albeit microscopic, 
value. 

There is actually a third possibility, that is in some sense half a way between the unbroken 
and broken symmetry case. We could allow for a time dependent bare mass, in such a way 
to simulate a sort of cooling down of the system. In order to do that, we could start with a 
unbroken symmetry bare potential (which fixes initial conditions naturally via the gap equa- 
tion) and then turn to a broken symmetry one after a short interval of time. This evolution 
is achieved by a proper interpolation in time of the two inequivalent parameterizations of 
the bare mass, eqs. ( |3T5| ) and ( |3.10| ). 

We looked for the infiuence this different choices could produce in the results and indeed 
they depend very little and only quantitatively from the choice of initial condition we make. 

As far as the numerical algorithm is concerned, we used a 4th order Runge-Kutta al- 
gorithm to solve the coupled differential equations (|A2|) , performing the computations in 
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boxes of linear size ranging from L = 207r to L = 4007r and verifying the conservation of 
the Wronskian to order 10~^. Typically, we have chosen values of corresponding to the 
UV cutoff A equal to small multiples of m for unbroken symmetry and of f a/A for broken 
symmetry. In fact, the dynamics is very weakly sensitive to the presence of the ultraviolet 
modes, once the proper subtractions are performed. This is because only the modes in- 
side the unstable (forbidden or spinodal) band grow exponentially fast, reaching soon non 
perturbative amplitudes (i.e. ~ A~^/^), while the modes lying outside the unstable band 
remains perturbative, contributing very little to the quantum back-reaction and weakly 
affecting the overall dynamics. The unique precaution to take is that the initial conditions 
be such that the unstable band lay well within the cutoff. 
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FIGURES 
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FIG. 1. Zero-mode amplitude evolution for different values of the size ^ = 20, 40, 60, 80, 100, 
for A = 0.1 and broken symmetry, with 4> = 0. 
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FIG. 2. Time evolution of the squared effective mass in broken symmetry, for ^ = 50 and 
A = 0.1. 
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FIG. 3. The quantum hack-reaction S, with the parameters as in Fig. 
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FIG. 4. Zero-mode amplitude evolution for different values of the renormalized coupling con- 
stant A = 0.01, 0.1, 1, for ^ = 100 and broken symmetry, with = 0. 
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FIG. 5. Detail of near t = tl for ^ = 40 (dotted line). The case ^ = 80 is plotted for 
comparison (solid line). 



3400 - 



3200 



2800 - 



2600 



2200 - 



2000 



100 



110 



120 



130 



140 



150 



160 



FIG. 6. Detail of ao near t = tl for ^ = 40 (dotted line). The case ^ = 80 is plotted for 
comparison (solid line). 
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FIG. 7. Next-to-zero mode (k = 2'!r/L) amplitude evolution for different values of the size 
■kz = 20, 40, 60, 80, 100, for A = 0.1 and broken symmetry, with = 0. 
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FIG. 8. Next-to-zero mode (k = 2it/L) amplitude evolution for different values of the renor- 
malized coupling constant A = 0.01,0.1, 1, for -kz = 100 and broken symmetry, with ^ = 0. 
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